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Abstract 

We consider elastic bodies in rigid rotation, both nonrelativistically 
and in special relativity. Assuming a body to be in its natural state in the 
absence of rotation, we prove the existence of solutions to the elastic field 
equations for small angular velocity. 
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1 Introduction 

The field of relativistic elasticity is still in its infancy. This work is part of a 
program where we set up the field equations and prove existence theorems for 
some of the most basic problems, both dynamical and time- independent. In 
the present paper we study the equilibrium of an elastic body in rigid uniform 
rotation. This is a time independent problem in the frame corotating with the 
body. Interestingly the nonrelativistic case of our result seems to be unknown, 
so we have to treat this also. 

We will be interested in equilibrium configurations of ideal elastic solids which 
are subject to the centrifugal force but otherwise free. The natural boundary 
condition, then, is that the so-called "normal traction", i.e. the components of 
the stress tensor normal to the surface of the body, be zero. The location of this 
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boundary of the region-in-space occupied by the elastic body can not be given 
freely, but is part of the sought-for solution. It is thus preferable to work in the 
material ("Lagrangian") representation where the maps describing configurations 
go from the material space - whose boundary is fixed - into physical space, with 
Neumann-type boundary conditions. We note in passing that it is the material 
representation which is used almost universally in standard nonrelativistic elastic- 
ity. Since the elasticity literature is hard to digest for workers with a background 
in relativity, we have made an effort to make this paper reasonably self-contained, 
by formulating the necessary concepts and computations in the framework of a 
Lagrangian field theory. 

Nonrelativistic elasticity in the time- independent case takes the form " 8,+H — 0" , 
where £ is a 2nd order partial differential operator acting on the configuration. 
This operator, usually called elasticity operator, is a quasilinear, elliptic opera- 
tor on Euclidean space, with coefficients depending on the elastic material. The 
force IF, usually called "load", depends on the problem at hand of course - in 
our case it is the centrifugal force. In special relativity this picture essentially 
survives with the complication that £ now becomes a partial differential operator 
(PDO) living on the space of trajectories of a "helical" Killing vector with the 
natural (curved) metric which this 3-space inherits from Minkowski space. For 
this reason - and for use in future work on elasticity in GR - we base our work 
here on a formulation of £ on the background of an arbitrary curved 3-space. 
In Sect. 2 we describe our setup for static elasticity, which is the curved-space 
generalization of the standard nonrelativistic theory for hyperelastic materials 
in the time independent case. The static elasticity operator is viewed as the 
Euler-Lagrange expression for a certain action principle. This is done merely for 
convenience, since the action principle facilitates the calculations necessary for 
moving back and forth between the spatial and the material picture. Wc also 
write down, in the material picture, the linearized elasticity operator at a "nat- 
ural state", i.e. at a solution of the field equations with zero stress (whence zero 
body force). When the physical space has Killing vectors, this operator, with the 
obvious choice of function spaces, is neither injective nor surjective. Rather its 
range has to satisfy certain integral constraints (often called "equilibration condi- 
tions" in the literature) on the force involving the Killing vectors and the natural 
configuration. There is a related fact concerning the full elasticity operator: any 
configuration and any force have to satisfy equilibration conditions, if the spatial 
metric has KiUing vectors. This lack of surjectivity of the elasticity operator 
generalizes the statement: forces acting on an otherwise free body at rest have 
to be such that the total force and the total torque be zero. This important fact 
is derived from the so-called balance laws in standard continuum mechanics. For 
a relativistic treatment see [11]. 

In Sect. 3 we derive the equations governing relativistic rotating elastic bodies. 
We do this by means of "dimensional reduction" of the general time dependent 
theory laid out in our previous work [2] , where this reduction is carried out w.r. to 
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the helical Killing vector corresponding to rigid rotation with angular frequency 
LO. The resulting action functional (rather: "energy functional") is in fact more 
general than the framework of Sect. 2: There automatically appears a force term 
in the form of a multiplicative function, namely the norm of the Killing vector 
which, in the case at hand, is essentially the centrifugal potential for frequency lo. 
As a complication, for purely relativistic reasons, the coefficients of the elasticity 
operator have a dependence on a;, which they inherit from the curved spatial 
metric, i.e. the natural metric arising by quoticnting the Minkowski metric by 
the action of the helical Killing vector. To avoid confusion we wish to stress that, 
under these circumstances, it is the full equation £ + 3" = 0, which is obtained as 
the Euler- Lagrange condition for the energy functional. 

Taking the formal limit c — oo of the field equations, one obtains the nonrela- 
tivistic equations. These of course have the form "flat elasticity operator + force 
=0", where the second term, i.e. the centrifugal force, is linear in uP' . In Sect. 3 
we solve these equations for small lo and configurations close to the natural one. 
An immediate application of the implicit function theorem is of course forbidden 
by the lack of surjectivity of the linearized elasticity operator in flat space. This 
is a well-known problem in elasticity, often resulting in bifurcation phenomena 
which have led to a lot of difficult work (see e.g. Sect. 7 of [9]). Our problem, 
luckily, turns out to be simpler. We first note that the equilibration conditions, 
when a; 7^ 0, require the configuration to be such that the center of mass in 
physical space lie on the rotation axis and that this rotation axis coincide with 
one of the principal axes of inertia: We ab initio impose these conditions on our 
allowed configurations (including of course the natural configuration, i.e. the un- 
deformed body) . The resulting space of configurations turns out to be a smooth 
manifold near the natural state (which it would not be if there were bifurca- 
tions). After a suitable projection of the equations to this manifold, reminiscent 
of bifurcation theory, the problem can then be solved using the standard im- 
plicit function theorem, provided the constitutive law (expressed in terms of the 
so-called "stored-energy function") satifies the condition of "uniform pointwise 
stability" , which is valid for standard elastic materials. 

In Sect. 4 we solve the relativistic problem. It has the form "quotient space 
elasticity operator + relativistic centrifugal force = 0" . For = the quotient 
space metric is Euclidean and the relativistic centrifugal force goes to zero like a;^ 
for small lo. If one now uses equilibration w. r. to "l/a;^x relativistic centrifugal 
force", a new situation seems to arise: for a; = we get the same conditions 
as before. But for cu 7^ we get none: the only Killing vectors of the spatial 
metric are now and c^a, and for those the equilibration conditions turn out to 
be automatically satisfied. Thus, in order to be able to use an implicit function 
argument at a; = 0, we resort to brute force: we split off the non-flat part of the 
elasticity operator, which vanishes like o;^, and view it as a contribution to the 
force. Now much the same goes through as in the nonrelativistic case. However, 
in order to be guaranteed that the set of equilibrated configurations is again a 
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manifold, we have to assume, in addition to the constitutive condition for the 
nonrelativistic case, that some characteristic velocity of the system, in essence 
some upper bound on the sound velocity, be sufficiently small compared to c. 
Being guaranteed, by the theorems of Sect. 4 and 5, that solutions exist for small a;, 
we finally, in Sect. 6, calculate these explicitly to linear order in tu^, for a material 
which is isotropic in its natural state and for an undeformed body which is an 
ellipsoid with rotational symmetry about the rotation axis. In the nonrelativistic 
limit our results agree with the ones found by [4], see [8]. 

2 The static elasticity operator 

The basic setup underlying time-independent situations in both nonrelativis- 
tic and relativistic elasticity is as follows: We consider maps between two 3- 
dimensional Riemannian manifolds given by / : {N,hij) i— > {Q,Vabc), with the 
manifold N describing physical space with smooth metric hij and Q a domain 
in (open, connected, bounded) with smooth boundary dQ (not necessarily 
connected) and Vabc a smooth volume form on Cl. The domain Q, called "body" 
or "material space", is to be t bought of as the collection of particles making up 
the elastic body prior to the action of any external forces, stresses, etc.^ Thus / 
is the "back-to-labels-map", its inverse ^ : Q N is called a configuration. In 
nonrelativistic elasticity [N, hy) is Euclidean space. (In the case studied here of 
a relativistic body rotating at angular frequency a;, the metric hij will be the one 
coming from the Minkowski metric on with the action of the helical Killing 
vector dt + ud^ quotiented out.) The maps / are required to be one-one and 
orientation-preserving, i.e. the function n, defined for each / by 

{dJ^){x){dJ^){x){duf){x)VABcU{^)) = n{x)e.,k{x), (1) 

is positive. The physical interpretation of / is that of the density of particle 
number. We are using coordinates X"^ on fl and coordinates on A^. The 
three form Eijk is the metric volume element in associated with and the 
three form Vabc{X) = V{X)eABC with V > and ei23 = 1 is the volume 
element on Q. (We think of V as having physical dimension [mass /volume].) 
Put differently the definition (Q) says that nh^ = Vdet{df), where h = det{hij). 
An elastic body will be specified by a constitutional law, as follows. There is given 
a scalar function of maps / called stored-energy function (of physical dimension 
[velocity]'^) w = w{f, df, x), smooth in all its arguments. Covariance of w under 

^It would be unnatural for our purposes to add further structure to the body at this stage - 
such as a flat metric, as is common in the literature. Once a choice of reference map ("state") 
has been made, there is of course a metric defined on by the push-forward of hij under this 
map. 



4 



spatial diffeomorphisms requires (see e.g. [2]) that w be of the form ^ 

w = w{H^^f), (2) 
smooth in its arguments, where 

= {dj^){d,nh^^. (3) 

By virtue of our assumptions H^^ is positive definite. It thus has an inverse Hab- 
The Cauchy stress tensor associated with w is defined as follows. Think of 
the function L = nw as a Lagrangian density, i.e. consider the action S[f] 



S = nwh^^ d^x. (4) 

J M 

Then is the Cauchy stress tensor 

d{nw) dw 
— Oij = I———. nwhij = ln——T (5) 

The mixed tensor is the same as the "canonical stress tensor" = nw6i^ — 
did^J^) ^^if^) corresponding to the energy functional (jH), which can also be 
written as 

We next define the elasticity operator £ by 

n£,[/] =D,a/, (7) 

where Di is the covariant derivative w.r. to hij. Clearly £j is a second-order 
quasilinear partial differential operator. Its relation with our variational principle, 
as can be checked in a straightforward fashion, is given by 

n£, = {dj^) &A, (8) 

where 8. a is the Euler-Lagrange expression corresponding to the energy functional 
S, i.e. 

, 1 /, 1 d(nw) \ dw 

Having written down the elasticity operator in the "spatial" (i.e. maps go from 
space into body) representation, we have to spell out the boundary conditions. 
For a free body these are the conditions of "vanishing surface traction", namely: 

cririjlf-^dQ) = 0, (10) 

^In standard nonrelativistic elasticity diffeomorphism invariance is replaced by requiring the 
validity of the so-called "principle of material frame indifference" (see e.g. [9]). 
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where Ui is an outward co-normal of the surface f~^{dQ). The awkward feature 
of the boundary condition Eq.()10|) of involving both the map / and its inverse 
$ = strongly suggests that one go over to the "material" representation which 
uses $ as the basic dependent variable. Note that such a change of representation 
should not be confused with a change of coordinates on either or Q. 
The field equations in the material representation are easiest to derive by again 
starting from the energy Eq.Q which now reads 

sm = [ wvd^x, (11) 

where the argument H"^^ in w should be interpreted as a function of ($, (9$). As 
such it can be viewed as the inverse of the pull-back metric 

Hab{X) = {dA^'){X){dB^^){X)h,,{^{X)) (12) 

of hij under $. In the elasticity literature Has is called the "(right) Cauchy- 
Green strain tensor" . Note that the material Lagrangian contains the dependent 
variable $ only through hij and has arbitrary dependence on the independent 
variable X, whereas the spatial Lagrangian contains the independent variable x 
only through hij and may depend on / in an arbitrary fashion. The material 
Euler-Lagrange expression, by virtue of Eq. (jSJ and the relation 

5^\f{x)) = -(9a$0(/(x))5/^(x), (13) 

is nothing but £j, expressed in terms of $ rather than /. Explicitly one finds 
that 

A calculation shows that Eq. p4|l can be written more concisely as 

= Vact,^, (15) 
where ai^ is the "first Piola stress tensor" given by 

and Va is given by 

(VAcr/)(X) = V-\X)dA{V{X)ai\X)) - V%[^{X))a^\X){dA^'){X) (17) 

with r*^ being the Christoffel symbols of the metric hij. Note that no metric on 
n is required, just the volume form Vabc The connection between the Piola and 
the Cauchy stress tensor is given by 

n{^{X))oi^{X) = ^^,{X)aimX)), (18) 
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where 

<^^,{X){dB^'){X) = 5^B. (19) 

or \l/'^j(X) = ((9j/^)($(X)). Thus the above derivation has in particular recov- 
ered a version of the so-called " Piola identity" , namely 



[n 



-^D,a,^){^{X)) = {VA<y^^){X) (20) 



The boundary conditions in the material picture take the form 

o"i = (Ti^nA\dn = 0, (21) 
where is an outward co- normal of dQ. From Eq. ()16|) it follows that 

From this we infer that the second-order terms of £i are given by 

S^a/ = [vl/^.ff^^i/cDa,^ + 2vl>^(,a,/ + 

+ 2H'''^H''^^'',ij'',LcEDF]{dAdB^n + Lo., (23) 

where ^ 

LabCD = I arrAR aTTrn 

(24) 



The problems considered in elastostatics are usually written as 

£i[$] + 3^,[$] =0 (25) 

The quantity is called "load". It will be convenient to use a more general 
terminology with respect to loads and elasticity operators. Namely, we call the 
(nonlinear) elasticity operator E the assignment, to any allowable map $, of the 
pair (£j, (jj), where (Tj is the function on dVL given by o"j = o-i^nA\dn- Similarly a 
load F is a pair of functions (?",, Tj) on Q x dVt, both of which may depend on $ 
(otherwise the load is called "dead"), and may do so in a nonlocal fashion. 
There is a condition on loads in order for Eq. ()25|) to have solutions, which will 
play an important role and which arises as follows: Let ^*(x) be a Killing vector 
of the metric hij on A^. Then, in order for a map $ to be solution of the extended 
version of Eq.(j2SI), i-e. E + F = the load F = (ff'j, Tj) has to satisfy 

{Co^)JidV+ [ (f o $) Ti = 0. (26) 
Q Jan 



The integrals in Eq. ()26|) ^ are with respect to the volume form Vabc on Q. The 
proof of Eq. ()26|) is a straightforward verification based on the Stokes theorem 



•^More precisely, the quantities at and should be interpreted as two-forms on dft with at 
corresponding to the pull-back to dfl of the two-form ai'^VABC- Similarly the first term in 
Eq. (|26|l should be interpreted as the integral over of the three-form arising by replacing 
with JiVABC- 
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and the following identity 



[dAiC ° + r}.(e o ^)idA^')]a,^ = 0, (27) 

which in turn follows from the Killing equation for ^ together with Eq. p8|) or 
Eq. ip^ (i.e. Eq.(j2Z|) is the material version of {D^^^)aij = 0). Loads satisfying 
Eq. (|26|) for all Killing vectors on $(f2) are said to be "equilibrated at Sim- 
ilarly, for given load L, a map $ is called equilibrated w.r. to L, if (I2(ij) is valid 
for all Killing vectors. 

We now introduce a reference configuration. This will simply be some given 
invertible, orientation-preserving map $ from fl to A^, obtained by restriction to 
of a smooth function defined in a neighbourhood of Q with smooth inverse on 

its image. We in addition require $~ to have n =P= const or, in other words, 
that 

V{X)=PH^X) (28) 

— — — — o 

with H = det(HAB), Hab being the pull-back under $ of hij. The constant P, 
usually called the mass density in the reference configuration, will in our conven- 
tion only appear in the expression Eq. (p?7j) for the elasticity tensor. The Eq. (fTH|l 
can now be written as 

e, = D^a/ (29) 

with Da the two-point covariant derivative (see [9]) referring to the metric Hab 
on VL and hij oji N ^. Given a reference configuration one can define, for any 
configuration the matrix Ti given by = HbcH^^. Much of elasticity 

theory concerns isotropic materials, which are described by stored-energy func- 
tions depending on H"^^ only via the principal invariants of Ti,. We will not need 
that assumption however. 

In this work we will always require that the reference configuration be unstressed, 
i.e. cjj^ = 0. By the identity unstressedness is equivalent to 

fd w{H<^^,X) \ _ 

gffAB ] \^=HX) - U- l^Uj 



The subscript $ = $(X) in Eq. (j3U|) is understood in the sense that $ and 9$ 
should be replaced by the values respectively of $ and 9$ at X. Clearly there 
holds 

= = 0. (31) 

We will need the linearization of the Piola tensor in the reference configuration. 
This is given by 

5a,^ = -2H^^^^,Lbcde SH^"" . (32) 

"^In contrast to [9] these metrics are not chosen independently, but are isometric under the 
reference configuration. 
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where Labcd = Labcb\>s>=^{x) ^iid SH^^ is the first-order hnear PDO acting on 
5$*(X) given by 

5/7^^ = Cs^H'^'', (33) 

and the Lie derivative in Eq. ()33|1 is taken with respect to the vector field 5$^(X) = 
^'^i{X)6^'{X) on Q. More explicitly Eq.(|22I) can be written as 

SH^^ = -2 ^(^(5$^) = -2¥^, D^^6¥, (34) 

since Da annihilates (see footnote ^•'). We remark that in the so-called "pre- 
stressed case" of linearization at a reference configuration which solves Eq. (l31|) 
without being stress-free there arises an expression more complicated than (jH^ 
for Sai"^ and hence for the linearized elasticity operator. 

One easily sees that 5H^^ vanishes if and only if is of the form (5$*(X) = 
^*(I>(X)) for ^* a Killing vector on ($(^2) C N,hij). The linearization of the 
elasticity operator E introduced after Eq.(j2SI) is the second-order operator 

6E:<^^ {DaScx,^, (5fx,V)|a^,) (35) 

with 5ai^ given be Eq.'s ()32p34|) . Clearly, functions {3^i,Ti) on x d^l, in order 
to lie in the range of 6E, have to be equilibrated at $. Now to the kernel. It is 
immediate that the kernel of 6E contains all elements 5$ which are "Killing" in 
the above sense or, equivalently, of the form (5$*(X) = {dA^^){X)ri^{X) , where 
r]^ a Killing vector on (0, Gab)- This is an expression, on the linearized level, of 
the following fact: given an unstressed state $, then, for any isometry F of A^, 
the map F o $ is also unstressed and hence a solution. 

We now suppose the stored-energy function w to be such that 6E is "uniformly 
pointwise stable". This means that, for X G fi, 

LABCDiX)M'^^M^^ > V M^^ = M(^^) ^ 0. (36) 

For example in the case of an isotropic homogenous material there are constants 
A and such that 

4po Labcd = ^HabHcd + '2fiHc{AHB)D, (37) 

and in that case uniform pointwise stability is equivalent to the conditions 

fi>0, 3A + 2/i > 0. (38) 

By a standard integration- by-parts argument, it follows from uniform pointwise 
stability that any element 5$ in the kernel of the map 6E is Killing. Thus 
the map 6E has a kernel of the dimension of that of the isometry group of 
($(f2),/ijj) and a cokernel of that same dimension. The latter fact can e.g. be 
inferred from the formal self-adjointness of SE, as follows: raising the index i 
in with h"^^ = k^^ o can be viewed as map from functions (5$*(X) 
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into themselves, which is formally self-adjoint w.r. to the inner product given 
by (5//, 5$) = J^hij{SiJ){6^^)d^V + j^^hij{5^x''){5^^)dO , using the symmetry 

LaBCB = LcDAB ^ ■ 

It follows from these observations, that, with an appropriate choice of function 
spaces, the operator 5L, viewed as a map from some complement of Killing ele- 
ments at $ to loads equilibrated at l», is an isomorphism ^. 



3 Rigidly rotating bodies 

In order to derive the equations governing a rigidly rotating elastic body wc start 
from the time- dependent theory as laid out in [2]. Back-to-label maps now go 
from a relativistic spacetime (M, g^j,) to (fi, Vabc)- The condition of invertibility 
in the time-independent case is replaced by the condition that there is a unique- 
up-to-sign timehke vector field with u^u'^g^i, — —1, so that 

u^id.f'') = (39) 

The quantity n, in the time-dependent setting, is defined by 

{d,f^){dj''){dxf)VABc = ne,,xpup. (40) 

The spacetime action for the elastic field is given by 

S= f ne{-g)'^d^x, (41) 



where e is a function of H"^^ = {d,j,f^){di,f^)g^^. We now suppose the maps 
/ are time-independent in the sense that is proportional to an everywhere 
timelike Killing vector field Furthermore suppose M — 'M} x N with N the 
quotient of M by the isometry group generated by The natural metric on 
N is given by hij = g.ij — goiOoj /Ooo in local coordinates (t,a;*) on x A^. Since 
gij _ f^ij ^ with h"^^ the inverse of hij, we have that the /'s, viewed as maps 
/ : {N,hij) ^ {Q,Vabc), and n bear the same relationship to each other 

as the quantities of the same name in the previous chapter. Since, furthermore, 
{~9)^ — {~9oo)^h2 we find that 

S ^ J ne{-goo)h^d^xdt, (42) 



^This formal selfadjointness runs under the name of the "Betti reciprocity theorem" in 
standard elasticity. 

^The necessary Fredholm theory is fairly standard in a Hilbert space (i.e. W^'^-)setting (see 
Ref. [13]). However in this paper we need W^'P with p > 3, and this can be found e.g. in [12] 
in the case when {N, hij) is flat. Luckily this is all we require for the present purposes. 
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and the reduced action reads 

S = I nt{H^^,f){-gooY^h^d'x. (43) 

We write e as 

e = c^ + w (44) 

and (700 as 

-^oo = c'e^. (45) 
The resulting field equations turn out to be equivalent to 

e-^Dj{e^ai^)-n{l + ^)DiU = Q (46) 

in the spatial picture and 



e-^VA(e^a/) - (1 + ^)D,U = (47) 



in the material picture. We now specialize to bodies in SRT which rigidly rotate 
at angular speed uj. Thus we take {M.g^y) to be {M.'^.rj^i,) with 



and ^^d^ to be 
There results 



rj^^dx^dx" = -c^dt^ + 6ij dx'dx^ (48) 
ed, = dt + ujd^. (49) 
e^ = l--^ (50) 



and 



,2 



hijix!'; u:)dx'dx^ = Sij dx'dx^ + ^ '~ 2 2 ^^^^^'^ ~ ^'^dx^Yi (51) 



UJ 

where = (x^)^ + (x^)^. We restrict ourselves to the region where r < -, i.e. 



inside the timelike cylinder on which gets null. Note that ^ enters Eq. ()46|1 . 
resp. Eq . fl47|) . explicitly, but ^ enters also implicitly through U and via the 
/lij-dependence both of the covariant derivative and of that of H^^ appearing in 
w, n and cTj-', resp. a^"^. 



4 The nonrelativistic Theorem 

Taking the formal limit of Eq. ()47|1 as c ^ 00 we find 

V~^dA{V(J,^)+u^{^\^\Q) = Q. (52) 
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Here cxj^ is given in terms of w by Eq-dTH)), with H^^ = \l>^j\l>^j S"-^ . Eq.^^ is 
the Euler-Lagrange equation for the action 

Snr = y K"^')' + {^^'fWVd'X (53) 

The action S^r is related to the action S in Eq. pHj) as follows. Take -S*, rewrite 
the integrand in material variables, discard a term proportional to which is 
a "null Lagrangian" i.e. does not contribute to the equations of motion and let 
c ^ oo: one obtains S^r- 
As usual the boundary condition is that 

or, = (aAA)lao = (54) 

Let $ be an unstressed reference configuration. It follows that $ is a solution 
of Eq. lfH^ and of the boundary condition Eq.(jnH). We can, and will, choose 
coordinates on VL so that $ is the identity map, i.e. = ^^{x'') = 6^iX^. 
From Eq. (PHj) we have that V = const in these coordinates. Thus, apart from an 
overall minus-sign, our field equation is of the form We are seeking solutions 
$ of Eq.'s (j52l54p for small values of u which coincide with $ for u = 0. We have 
the conditions (|26|) for each element of the Lie algebra of the Euclidean group. 
These six conditions, in the presence of fl54|) . amount to the statement that, for 
any configuration $ solving the field equations and the boundary condition, the 
total centrifugal force and the total centrifugal torque be zero. The former of 
these conditions, namely 

co^ [ $^3^ = u;^ [ <l>^d'X = (55) 

states that the center of mass has to be on the axis of rotation. The vanishing of 
the total torque reads 

[ <^^<^^£X = uj'^ ! <l>3$2rf3x = 0, (56) 
Jn Jn 

and this means that the axis of rotation is an eigendirection of 

Qij = j _ $^$J)rf3x, (57) 

Jn 

apart from a factor j^pd'^X the tensor of inertia. Hence rigid rotation is only 
possible through a principal axis of inertia of the rotating object. 
The above conditions depend on the configuration $ we want to determine. If 
one linearizes at the stress free configuration a kernel and range of the linearized 
elasticity operator appears, as explained in Sect. 2, and one can not use the im- 
plicit function theorem directly. In the known existence theorems this problem 
was dealt with as follows: 
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Stoppelli [10], who gave the first existence theorem for dead loads, in a first step 
projects the equation on the range of the hnearized operator. In a second step he 
shows that one can use the invariance the hnearized operator under motions to 
construct, from the "projected solution" , a solution of the original equation. This 
is similar to the method of Liapunov-Schmidt reduction in bifurcation theory. By 
a similar technique, Valent [12] solves various problems including ones with life 
loads. 

We propose here yet another method which rests on the fact that our life load is a 
differentiable function of the configuration. Therefore we can construct a priori a 
manifold of configurations which are equilibrated for the centrifugal force. Using 
only these configurations the implicit function theorem directly gives the solution 
for small cu, provided our stress free initial configuration $ is equilibrated for the 
centrifugal force and has three different moments of inertia. 

We now fix the spaces we will work in to be the standard ones in static elas- 
ticity problems. We assume that the boundary of our body Q is at least C^, 
i.e. there are no corners. We take the components of the configurations in 
W'^'^{fl,M.^) for p > 3. Then $ is in C^(f2) and the boundary traction cTj'^n^ is 
in W^-^/P'P{dn,R^). 

We can choose a neighbourhood C of $ in W'^'P{Q, M^) small enough so that each 
$ G C has a C^-inverse (see p. 224 of [5]). The elasticity operator E of Sect. 2, 
given by 

E : $ G e ^ [dACTi^, (Ti) G W^'P{Q, R^) X W^-^/P'P{dQ, M^) (58) 

is differentiable. (For the first factor in see Appendix A, for the second 

factor, see Remark (6.5) on p. 78 of [12].) We call 

L = w°'P{n, w^) X w^-^/P'P{dn, w^) (59) 

the "load space" and its elements (3'j,rj). The load map F : C — > L is given 
by = uj'^Zi together with = 0. The form of the centrifugal force Zi = 
— ($^, 0) and the principle of material frame indifference imply (see footnote 
^^) that, given any solution ($^,$^,$^) oi E + F = 0, ($\$^,$^ + c) and 
(<|)'i^<|)'2^$3^ where ($'\ $'^) are a rotation of ($\ <l>^) in the (x\ x^)-plane, are 
also solutions. We fix this freedom by imposing the conditions 

6° := {$1$ G e, $=^(0) = 0, (92$^)(0) = (9i<l>2)(0)}, (60) 

where we have also assumed that Q has been chosen such that G ^2. Clearly C*^ 
is a C^-(in fact: analytic) submanifold of C. Finally we want to restrict ourselves 
to configurations, which are equilibrated w.r.t. the centrifugal force, i.e. we define 

6° := {$1$ G e°, [ <i>^d^X = [ (l>^d^X = [ <^^<i>^d^X = [ <^^<^^d^X = 0}. 
Jn Jn Jn Jn 

(61) 
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Suppose furthermore that $ G C^;- Exphcitly this means that 

/ X^d^X = I X^d^X = I X^X^d^X = [ X^X^d^X = (62) 
Jn Jn Jn Jn 

Note that the last two equations in ()62p mean that the rotation axis concides 
with some fixed but arbitrarily chosen principal axes in the given reference con- 
figuration. Then we want to show that is a submanifold of finite codimension 

in e°. 

The map H, which sends a configuration to the values of the four integrals in 
Eg . (1611) is a differentiable map from C° to M^. Hence C C° has codimension 
4 near $ if the four differential 1-forms defined by the four conditions above are 
linearly independent at $. These 1-forms are 

A(5$) = [ 5<^^d^X , 5(5$) = [ 5<^^d^X (63) 
Jn Jn 

and 

C(5$) = [ {X^5^^ + X^5<^^)d^X , D{5<^) = [ (XW + XW)rf^X (64) 
Jn Jn 

Suppose a linear combination exists, i.e. 

aA{S^) + bB{6^) + cC(5$) + dD{6^) = (65) 

for constants {a,b,c,d) E R.^, which vanishes for all 5$* G PF^'^(n,R^) in the 
tangent space of at $, i.e. satisfying 6^^{0) = 0,{di6^^){0) = (925$^)(0). 
Choosing 5$^ = 5$^ = we obtain 

[ {cX^ + dX^)6<!>^d^X = 0. (66) 
Jn 

Choosing 5$^ = X^ and successively = and using the Schwarz inequality, 
we infer c = d = 0. Similarly we obtain a = b = from considering 5$^ = (X^)^ 
and 5$^ = (X^)^. Maximal codimension at $ implies maximal codimension 
nearby. Therefore the level sets of H are closed submanifolds with the property 
that the tangent space splits into an infinite dimensional closed subspace tangent 
to H = const and some finite dimensional subspace (see e.g. Theorem 3.5.4 
of [1]). We note that the above proof works for any $ which is equilibrated. No 
condition like "no axis of equilibrium" (see [12]) appears at this point. We collect 
the above rsults in the 

Lemma 1: The configurations in C^; form a differentiable closed submanifold of 
codimension 4 in C° for $ sufficiently close to the identity. 

We can say more about possible finite dimensional complements of C^,. Let us 
assume the equilibrium configuration {Q, $) to be such that all three moments 
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inertia are different, hence the principal axes of inertia are uniquely determined. 
This implies that any small rotation around the or - axis and translation 
in x^, - directions will destroy equilibration! Therefore the corresponding 4 
Killing vectors span a complement of at $. 

We now turn to the nonlinear map 

S-.elxR^L (67) 

defined by the l.h. side of Eq.(jH21), i.e. 

g[$;o;] = (£4$] a,\^]) (68) 

We have that Si^*, 0] = 0. We want to solve the equation = for small u 

by the implicit function theorem. The map S is differentiable (see Appendix A). 
The derivative of Q with respect to $ at $ and c<j = is the linearized elasticity 
operator together with the linearized normal traction on the boundary, i.e. 

^[*;^]l(*;0)= (^[<^>]|*,^[<^]|#) =5^:7^#(e°z)-i^. (69) 

We now recall the discussion of Sect. 2. If we consider the elasticity operator 6E 
as a map C ^ L we have KeT{6E)= Killing vectors and Range(5£^) = L^, the 
loads equilibrated at $. Hence on the tangent space of 6° the kernel is trivial 
as we showed above. To deal with the range we proceed as follows: Choose a 
(6-dimensional) complement 5*6 to Lj, which defines a unique linear projection 

P:L = L^®Se^L^ (70) 

As is done in bifurcation theory we can solve the "projected equation" i.e. 

P o (£,[$]+ ^2^,;[$]) =0, Poa4<l>]=0 (71) 

by the implicit function theorem for small u ^ because the derivative of P o E, 
namely PoSE, is an isomorphism at 5>. Denote this configuration by Finally 
we show that already solves our problem: fix any configuration $ G and 
consider the codimension-6-linear-subspace L$ C L of all loads equilibrated at 
$. These are all (3",, r^) satisfying 

[ 3'i£X+ [ TidO = , [ A 3^)i d^X + [ ($ A r)i rfO = (72) 
Jn Jdn Jn Jan 

The projection P, when restricted to L<j,, defines an isomorphism P^ : Lq, ^ 
(see [7]) because, for configurations $ near $, the complement Sq of is still a 
complement of L$. In particular L<j, and just intersect at the origin. For all 
$ G we have 

(£,[$], (T,[$]) gL$ (73) 
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Hence this holds for Since $a; £ there holds 

(Z,[<l>J,0)eL$„ (74) 

by construction of C^; • Hence we obtain 

+a;2Z,,[<l>J,aa<f J) G (75) 

Our candidate solution $^ satisfies Po(£j[$^]+co'^2.j[$t<;]) = and Pocrj[$^] = 0. 
Consequently 

£,[<I>J + cu2Z,;[$^] = 0, = 0. (76) 

Hence $^ is a solution, and we have proved the 

Theorem 1: Let the domain Q and the map $ be such that all three moments of 
inertia are different and that one of the principal axes concide with the rotation 
axis. Let the stored-energy function w be uniformly pointwise stable, i.e. satisfy 
the inequality (|36|) . Then there is, in a neighbourhood of ($, a; = 0) G C° x R, a 
unique element in C^, solving the equation together with the boundary 
condition (j^H). 

We remark that the condition on the natural configuration, in Theorem 1, of 
having three different moments of inertia /i, /2, /a, is not necessary for the theorem 
to go through. In Sect. 7 we consider a natural configuration which, by virtue of 
its axisymmetry with respect to the rotation axis, has Ii = I2 I3, with the 
eigenvectors for Ji and I2 orthogonal to the rotation axis. Then Theorem 1 
remains true, with the proviso that the rotation axis is chosen to coincide with 
the I3- axis, since the only additional freedom of performing continuous motions 
is that under d^p and d^, and that has already been frozen out in C^. 
Consider finally the case of a spherical top, i.e. where all eigenvalues of the 
tensor of inertia coincide. An example is given by taking Q to have the geometry 
of a cube (w.r. to the the metric Hab)- Our above proof works also in this 
case: We now have a 2-dimensional intersection of the kernel of the linearized 
elasticity operator with the tangent space of C^, at We can solve by the implicit 
function theorem if we fix an element of the kernel. The cube is equilibrated for 
the centrifugal force for any rotation axis which goes through its center of mass. 
Hence there is a 2-paramter family of possibilities for $. For the cube these $'s 
will determine physically distinct non-linear solutions. On the other hand, for a 
sphere, the trivial spherical top, different choices of $ lead to nonlinear solutions 
which are just rotations of each other. In general, the interplay between the 
symmetry of the tensor of inertia and the symmetry of the domain Q determines 
which of the solutions we obtain by selecting an element of the kernel are different 
Remark: If the stored-energy function is analytic, we obtain maps between the 
function spaces which are also analytic. The analytic implicit function theorem 
implies that the family $^ is analytic in a;^, i.e. we have a converging Signorini 
expansion. Furthermore elliptic regularity in that case implies analyticity of the 
solutions in Q. 
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5 Non relativistic, self gravitating, rotating 



Combining the result of the last section with our paper on a static self gravi- 
tating body [3], it is straightforward to obtain an existence theorem for a slowly 
rotating, weakly gravitating body. Namely, if we add the gravitational force to 
the centrifugal force, the equation to be solved takes the form 

The important point, then, about self-gravity is that the last term in Eq. (f77j) is 
automatically equilibrated for all configurations. One can then basically proceed 
as in the last section ''. 



6 Relativistic case 

A specific model is again characterized by a choice of stored energy w = w{H^^, X) . 
The main complication is that H^^ now refers to a curved spatial metric which 
depends on u. In moving back and forth between the relativistic and the non- 
relativistic theory, it is natural to require that the stored-energy functions, for a 
given material, be the same. It follows that 



fdw{H^^,X) 



\(<S>=<S>(X);uj=0) 



(7J 



and that 

^ABCD - 1 dJjABQJjCD J I (*=*{X);<^=0) ['^) 

be pointwise stable in the sense of Eq. p6|) . As opposed to the nonrelativistic case 
the equations are not invariant under adding a constant to w. We thus assume 
that 

w(i/^^,X)|(^=3,(^);,=o) =0 (80) 

As in the previous section we assume coordinates X^ to be chosen such that $ 
is the identity map. In the relativistic case this requires that [X^Y + (X^)^ < 
for all (X\X2,X3) e n. 

We will treat the relativistic case by splitting off the nonrelativistic elasticity 
operator in Eq. ()47p and putting all remaining terms into the load. This we do as 
follows: We can write 

e^a/ =a,^ + u^^, (81) 



^We take this opportunity to point that, on p. Ill of [3], we required the stored-energy 
function to satify the so-called Legendrc-Hadamard condition, whereas we should have required 



the stronger condition H36|l employed here, in order for the theorem stated there to be true. 
Furthermore the estimate a 
M is the mass of the body. 



\T\L 

Furthermore the estimate after Eq.(2.1) of that paper should be replaced by G < ^jp- where 
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where = 9$; lu) and =o"i^U=o, and analogously 

{e-^a^^nA)\an=^i+^% (82) 
for the surface traction. The boundary conditions now take the form 

-00% (83) 

By virtue of Eq. (f7S|) we have that 

^^^m = 0. (84) 
The field equations (j47j) can be written as 

dAa,^ + u;\'^i + Zi) = 0, (85) 
where ^ 

= 9a?/ - a f 1 - ^1 ' (5A$^)rt. a/ (86) 

and 



where it is understood that the functions r and V'l- are evaluated at the points 
$(X) G A^. Note that the second term in Eq.(|HU|) is regular also at = 0, due 
to Eq.(jnH). Furthermore, using and (jHDI), there holds 

Za<l>;cu]|($,o) = (X\X2,0). (88) 

We want to solve Eq. ()85|) by viewing the first term as the (unperturbed) elasticity 

o o 

operator £j. (Note that the fiat-space operator £j, though "unperturbed", is still 
nonlinear.) The remaining terms in Eq.()85p form the load, i.e. 

= u\Z, + ^i) (89) 

together with 

Ti = uJ%. (90) 

Thus the equilibration conditions Eq. lj^Bj) . dividing by cu^ and turning the surface 
integral to a volume integral, take the form 



= - / dA^i'o ^)a,^d^X + 

+ / {Co^)Z^[^;uj]d!'X. (91) 
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The vectors in Eq.()9ip run through the Euchdean Kilhng vectors on A^. We 
now claim that, as in the nonrelativistic case, 2 of these 6 conditions are identities, 
namely, if ^* is either or d^f,. To see this, recall that the conditions (j9T| are 
equivalent to the relations 

/ C [—oDAai^ + Zi] £X = (92) 
Jn ^ 

By virtue of the axial symmetry of r^, the second term on the right in ()92|1 gives 
zero when ^ is c^s or dcf,. The first term, using the boundary conditions and that 
^3 and are Killing vectors also of the c<j-dependent curved metric /ij-,-, gives also 
zero for 7^ 0, whence for = 0, by continuity. This proves the above claim. 
We next have to look at the different terms in the integrand of Eq. ()9ip at u; = 0. 
We have that 

,2 



W IT 

Z.U=o = (l + ^)9j - ) . (93) 



Furthermore we find that 
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^r,.cr,, U=o = Ok , (94) 



where the quantities Af^ are linear functions of $(X) with constant coefficients 
and, using Eq.'s (f^OI^IIHI jl and Eq.dSl, 

?/U=o = ^/ - I k^^ a,^ kBC -| iscDE, (95) 

where K"^^ = ^ H^^\^=o = '^'^i'^^ jK."-^ and k*-' are quadratic functions of $(X) 

o o 

with constant coefficients, Hab= Hab\uj=o and Labcd= Labcd\uj=o- Note that 
all quantities with superscript ° depend on d^{X), but not on ^{X). We next 
evaluate the equilibration conditions in the reference configuration. The first 
term in Eq.(|91|) gives no contribution at $ = $, since the first two terms in (j95|) 
vanish and the third term contributes zero, due to the symmetries of Labcd and 
the Killing equation for ^. The second term in (j^Tj) is also zero in the reference 
configuration. The third term, finally, is identical with its nonrelativistic value, by 
Eq. ljHUj) . It follows that the requirements in the previous section on the reference 
configuration can remain unchanged. It remains to compute the derivative at 
$ = $ of the function given by the four functions resulting by inserting 

into the r.h. side of (j9ip the Killing vectors ^ = di,^ = 82-,^, = x^di — x^d^jC, = 
x^d2 — x^ds on A^. The derivative at $ = $ of the last term in ()9H). using ()93|) . 
the stressfreeness of $ and ()80|). is the same as that in the nonrelativistic case. 
The explicit form of the remaining terms does not matter except that they are 
linear, with coefficients some given functions of X, in the quantities -^Labcd{.X). 
The quantities -^Labcdef{X), where 

^'"^ ' ' (96) 



ABCDEF - 1 gjjAB QffCD QffEF ) I [* =1'(X); <^=0] 
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appear in the derivative of Eq.(|U3j). but do not contribute to DH at $ = $, again 
using the symmetry and the flat-space Kilhng equation for ^. Now recall the 
discussion of Sect. 4. Evaluating the nonrelativistic DH\^ on some test functions 
5$, we arrived at certain positive expressions. These can bounded from below 
by a (dimensionful) quantity, say (3, which solely depends on the geometry of 
the domain VL. Similarly, the contribution to DH\q, of the relativistic terms just 
discussed can be bounded from above by a geometrical quantity, say 7 times 
^|L|, where is some upper bound for the components of Labcd- It follows 
that there is a dimensionless number a, which depends only on the geometry of 
so that Dif |[$^,|. ^^0] is nonzero provided that 

^ < a = ^ (97) 

We now follow the pattern of the discussion in Sect. 4 as much as possible. We 
assume the domain VL to lie strictly inside the cylinder {X^Y + {X'^Y = %• We 
choose a neighbourhood '^D of the identity in C, small enough so that C iV^, 
for all < cij < e. Here is the subset in with r < -. We then restrict 
^2) to "D", by imposing the conditions Eq. ljUUj) . The elasticity operator E is still 

the nonrelativistic one, namely $ 1-^ (£j =£»= — SaC"?/^, c^i). The load map F, 
consisting previously of J'j = ci;^($^, 0) together with = 0, is replaced by 
J. = cu^Z', = Lu'^C^i + Z,), with Z,^ according to Eq.'s (I8eil87jl . together with 
Tj = ciJ^CTj. The important difference is that Z'^ depends on uj and on (9$, (9(9$) 
with (99$ appearing only linearly as required by the second result in Appendix 
A. Let us denote by ^D^./ x (— e, e) the set of configurations in ^D'^ and values 
UJ G (— e, e) satisfying the 4 relativistic equilibration conditions (i.e. the zero-level 
set of the function H described above, which involves Z'^). Our above discussion, 
together with the inverse function theorem, shows the following 

Lemma 2: Suppose the inequality ^7\i is vahd. Then the set "Ti^^, x (— e,e), 
for $ sufficiently close to the identity and u sufficiently close to zero, is a C^- 
submanifold of codimension 4 in the Banach space x R. 

We now consider the equation Po[E+F) = on "^T)^, x (— e, e) with the projection 
map P defined exactly as before. The remainder of the argument is completely 
analogous to Sect. 4, and we obtain the 

Theorem 2: Let the stored-energy function w satisfy Eq. (j8U|) and Eq.(|78p. and 
let the constants Labcd defined by (f?^ satisfy the inequahties (jHUj) and (jHIj). 
Suppose, finally, {Q, $) to be such that the three principal axes are different and 
that one of them coincides with the rotation axis. Then there is, in a neigh- 
bourhood of ($, 0) e "D% X (— e, e), a unique element ($^,ti;) which solves the 
equations (jH^ together with the boundary condition (jH^ . 
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7 Linearized solutions 

In this section we want to present an explicit solution which is linearized in lo^ . 
We take as the background solution a stress free ellipsoid of the form 

+ + e^Z^ = (98) 

The material is assumed to be isotropic in its natural state. Thus we have that 
Labcd defined by Eq. (f7^ satisfies (see Eq. (IHZI)) 

4po Labcd = ^^ab^cd + '^fJ'Sc{ASB)D, (99) 

and the inequalities PHj). namely 

/i>0, 3A + 2/i>0. (100) 

Consider the family of solution of (82) determined by the implicit function theo- 
rem and parametrized by u"^. The linearization in u"^ satisfies the equation 

a^^aZ + y. + z, = 0, (101) 

where, again, the bar means evaluation at $ = $ and u = 0. 
We write 

= ^"^-ic-o (102) 
and (5a/ is determined from dSSESEl)- We obtain, using (inni),(IZSl)and that 

6a,A= -H^AS^^ + d,6^A + XS.AdtS^'']. (103) 
Po 

Note that, by the convention to view $ as the identity map, = and 

Hab = Sab- Eq. ()l()3|l leads to the standard nonrelativistic operator of linearized 
elasticity, i.e. 

Oa S^, ^ = -[/iA<5$, + (/i + X)dAS^'']. (104) 
Po 

For the remaining terms in the equation we obtain from (j86l87|) 

^^ = ^{\-li){X\X\0) (105) 

Zi = {X\X\0) (106) 
The boundary conditions at (X^)'^ + (X^)^ + e'^{X^Y = are 

(5a/ + ?/)nA|aQ = (107) 
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where we can take n^rfX^ = X-^dX^ + X'^dX'^ + e^X^dX^ and has to be 

evaluated at $ = $ and u = using Eq. (jU3j) . 

We can replace A by the Poisson number a of elasticity defined by 

A = ^ (108) 

which satisfies 

- 1 < (7 < - (109) 

Then /i drops out from the boundary conditions and in the linearized equations 
/X appears only as — . As the "force is equilibrated", the linearization of the 
boundary value problem has a unique solution satisfying 5$^(0) = 0, (925$^)(0) = 
(9i5$^)(0) = (see Eq. ()60|) ). For the non relativistic case a solution of the above 
equations can be found in [8] going back to [4] . We make to following ansatz for 
6¥: 



6^' = X'[ai + a^iiX'f + (X^) + asiX^'f] (110) 

<5$2 = X2[ai + a2iiX'Y + (X^) + asiX'f] (111) 

5<l>3 = X3[a4 + a,i{X'Y + (X^)^) + ae{X^f] (112) 

The rational behind this ansatz is that we have a linear PDE-problem with con- 
stant coefficients, a polynomial (in factdinear) right-hand side and boundary 
conditions on an algebraic surface. Thus the solution should also be polynomial. 
In fact, inserting the ansatz Eq. pi0lllllll2|) into the equations and the bound- 
ary conditions we obtain by a lengthy calculation a linear inhomogeneous system 
Qa = C for a = {ai,..aQ) (note that 2 of these equations come from the field 
equation and 4 conditions come from the boundary conditions) . Using Maple we 
find 



Q 



2Xe\AXR'^e'^,2fiR\{2n + X)e'^,{4fi + X)R^e'^,0 
0, -4A£^ (-2/i + 2A)£2, 0, (-2/i- A)£^-2/ie2^ (6/i + 3A)e2 

0, 0, 4A + 4/i, 0, 4/i, 12/i + 6A 
0,16/i + 8A,2/i,0,2A + 2/i,0 
2A + 2/i,4Ai?2 + 6/ii?2^0,A,Ai?2,0 
0, (-6/2 -4 A) 4/i + 2 A, 0, (-A + 2 /i) 3 A 



with determinant 



det{Q) = -48 (2 /X + 3 A) (2 /i + A) (6 £ V + 64 /i + 11 A + 64 A + 20 A e^) 



and 



2c2 ' 2c2 ' 
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Here are some observations concerning this linear system: 

1. ) For positive /i, A obeying ()100p . we can solve the linear system. 

2. ) The velocity of light, c, appears only on the right-hand side of the linear 
system and we can write the solution as the non-relativistic solution plus a rel- 
ativistic correction term proportional to ^. For c ^ oo the vector C greatly 
simplifies. 

3. ) The case A = 0, i.e. a = 2X^2/!' much simpler. In particular, the only 
relativistic correction is a change of the "effective density", — po —Po + c^^ifJ' — 
A). 

The solution can be given in closed form. We begin with the simplest case, i.e. 
e = 1: the deformation of a sphere. We find that 

_2 {-3 + 2a + 3a')R^ po 

5 (a-l)(5a + 7)(a + l) ^ ^ ^ 

1 (-4 + 3a + 5a^) po 
"^-"10 (5a + 7)(a-l)7 ^ ' 

1 (-9 + 8a + 5a^) po 
"^=^-"10 (5a + 7)(a-l)7 ^ ' 

1 (-3-18a + 3a^ + 10(T^)i?^ po 
"""" 10 (a-l)(5a + 7)(a + l) ^ ^ ^ 

„ _ 1 - 3) Po 

"'-5(5a + 7) {a-l)J 

ae = -- , + (118) 
^ 10 (5 a + 7) ((7- 1) /i ^ ^ 

The change of a point X is given by uj'^6^^{X). On the equator and north pole 
one finds: 

5$^(0, R, 0) = Ria, + a,i?^) = -1 (119) 

10 (a + 1) (5 (J + 7) /i 

**»(0, 0. fi) = fl(a, + a,R^) = -1 I^^^-tl^Llll^^* . (120) 

5 (cr + 1) (5 cr + 7) jj 



Note that the north pole can move outward when — 1 < cr < — ~^ , but for 
physical materials a will be positive. Next we give the relativistic corrections 
which we denote by bi . They are independent of yu: 

3 i?2(5a2 + 3cT-4) 
bi = -tt: 



b2 = ^ 



10 (5 a + 7) (a - l)c2 

1 5^2 + 11(7-4 



10 (5 a + 7) ((T-l)c2 
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64 = -tt: 



3 10 - 3 (J - 3 
10 (5 (7 + 7) (cr - 1) c2 

3 (3(7 + 1) 



10 (5 (7 + 7) ((7 - 1) c2 

3 5(72-7(7-2 

10 (5 (7 + 7) ((7 - 1) c2 

1 1+ (7 +10(72 

^ ~ 10 (5(7 + 7)((7- l)c2 

The change of a point on the equator and the north pole are: 

2 (15(7=^ + 6(72 - 17(7 + 2) 
5$ = 

5 (-l + 2(7)((7+l)(5(7 + 7)c2 

3 _ _ 1 (30(7^ + 47(72 - 4(7- l)i?3 
~ ~5 (-l + 2(7)((7+l)(5(7 + 7)c2 

To show the effect of the eUipticity, we give just oi : 

PqR^ (6 + 3 £2 + 3 £^ + 30 (72 + a2 ^ ^4 ^3 _ 27 ^2 ^2 ^ 23 - 5 a - 28 a + a e^) 

/X ((7 - 1 ) ((7 + 1 ) ( 10 £4 + 40 (72 £2 _ 8 0- ^2 + 48 ^ + 5 ^4 _ 15 _ 1 1 £4 _ 8 £2) 



The hmits e = 00 exist. They approximate a very flat elhpsoid. The formulas 
are comparable to the nonrelativistic case. 
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A Some Functional Analysis 

We collect some differentiability statements which can be easily extracted from 
Ref. [12]. First consider maps / of the form 

/ : r3 X ^9 ^ ^3 

Let $ G W^'P{n,R^), p > 3. Define the "Nemitsky" operator / : W^'P{n,R^) 
W^'P{n,R^) given by 

M^){X):=UHX),d^X)) (2) 
Then / is if / is and / is C"^ if / is C"^. 
Secondly, consider maps g of the form 

^ : X ^ X (3) 
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Then the quasihnear operator g 

gm{X) := gtl{^{X),d^{X))d,di^\ (4) 

viewed as a map g : W'^^^{Vt,R^) W^°'f(Q,R3), is if / is and if / is 

An elementary statement used in the body of the paper is the following: If $ is in 
LP with p > 1, then the map fi : W — > M sending $ into $ d?X is continuous, 
thus analytic by linearity. 
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